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Abstract 

We have examined the importance of momentum dependent induced nucleon 
currents such as weak-magnetism and pseudoscalar couphngs to the amph- 
tude of neutrinoless double beta decay (Oi//3/3-decay) in the mechanisms of 
light and heavy Majorana neutrino as well as in that of Majoron emission. 
Such effects are expected to occur in all nuclear models in the direction of 
reducing the light neutrino matrix elements by about 30%. To test this we 
have performed a calculation of the nuclear matrix elements of the experi- 
mentally interesting nuclei A = 76, 82, 96, 100, 116, 128, 130, 136 and 150 
within the proton-neutron renormalized Quasiparticle Random Phase Ap- 
proximation (pn-RQRPA). We have found that indeed such corrections vary 
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somewhat from nucleus to nucleus, but in all cases they are greater than 25%. 
In the case of heavy neutrino the effect is much larger (a factor of 3). Com- 
bining out results with the best presently available experimental limits on 
the half-life of the Oi^/3/3-decay we have extracted new limits on the effective 
neutrino mass (light and heavy) and the effective Majoron coupling constant. 



PACS numbers:23.40.Hc,21.60.Jz,27.50.+e,27.60.+j 



I. INTRODUCTION 

The neutrinoless double beta decay (Oz//5/5-decay) is expected to occur if lepton number 
conservation is not an exact symmetry of nature. It is thus forbidden in the Standard 
Model (SM) of electroweak interaction. The recent Kamiokande results have given evidence 
that the neutrinos are massive particles and one has to go beyond the SM. To further 
understand neutrinos, we must know whether they are Dirac or Majorana particles, an issue 
which only double beta decay can decide. The Oi//?/?-decay can be detectable only if the 
ordinary beta decay is forbidden or suppressed and the neutrino is a Majorana particle (i.e. 
identical to its own antiparticle) with non-zero mass The study of the Oz//5/5-decay 

is stimulated by the development of grand unified theories (GUT's) and supersymmetric 
models (SUSY) representing extensions of the SU{2)l ® U{1) SM. The GUT's and SUSY 
offer a variety of mechanisms which allow the Oi//?/?-decay to occur |Q. The best known 
possibility is via the exchange of a Majorana neutrino between the two decaying neutrons 
but increased attention is paid to more exotic processes, like the supersymmetric 
R-parity violating mechanisms of Oz//5/5-decay Recent review articles |TB|JT^ give a 



detailed account of the latest developments in this field. 

In this contribution we shall discuss the role of induced currents such as weak-magnetism 
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and pseudoscalar coupling in the calculation of the Oz//3/3-decay amplitude, which enters the 
neutrino mass as well as the Majoron emission mechanisms. So far only the axial-vector and 
the vector parts have been considered systematically and in great detail p|-^j6| JT^JI^ . 



The weak-magnetism and nucleon recoil terms have been considered in the extraction of 
the neutrino mass independent parameters associated with right-handed current mechanisms 
of Oi//3/3-decay and were found to be very important in the case of the rj parameter []6|,[l5 



This is understood since the leading contribution in this mechanism is proportional to the 
lepton momenta. In the two-neutrino double beta decay (2z//5/5-decay) the weak-magnetism 
term |jl6| resulted in a renormalization of the Gamow- Teller matrix element independent of 



the nuclear model and lead to a reduction of only 10 percent to the half-life of medium heavy 
nuclei. In the case of the light neutrino mass mechanism of Oz//?/3-decay there has been one 



attempt to include such effects resulting |17] from the recoil term and their contribution 



has been found very small. The weak-magnetism and induced pseudoscalar terms have been 
considered in connection with the heavy neutrino mass exchange mechanism in Ref. and 
their importance have been manifested for Oi//5/?-decay of ^^Ca. 

To our knowledge the induced pseudoscalar term, which is equivalent to a modification 
of the axial current due to PCAC, has been ignored in all calculations studying the light 
Majorana neutrino mass mechanism even though it provides a contribution, which is in 
fact greater than, the included in all calculations, vector contribution. This has perhaps 
happened because in the charged- current weak processes the current-current interaction, 
under the assumption of zero neutrino mass, leads to terms which except the vector and 



axial- vector parts |]T9[ are proportional to the lepton mass squared hence, i.e. they are small. 

The induced pseudoscalar term, however, is a real function of the Lorenz scalar there- 
fore there is reason to expect it to be important. In fact we find that such corrections are 
of order {q)'^/{{qf + m^)^ i.e. they are important if the average momentum < g > of the 
exchanged neutrino can not be neglected in front of the pion mass. In the case of a light 
intermediate neutrino the mean nucleon-nucleon separation is about 2 fm which implies that 
the average momentum < g > is about 100 MeV. This leads to corrections of about 30%. 



In the case of a heavy neutrino exchange the mean internucleon distance is considerably 
smaller and the average momentum < g > is supposed to be considerably larger. 

We should mention that in the R-parity violating SUSY mechanism of Oi//3/5-decay ||Tl 



one has scalar, pseudoscalar and tensor couplings at the quark level, which, of course, induce 
analogous couplings at the nucleon level. 

The correct nucleon current is important in any calculation of the nuclear matrix ele- 
ments, which must be computed precisely in order to obtain quantitative answers for the 
lepton number violating parameters from the results of 0z//3/3-decay experiments. 

The goal of the present paper is to obtain reliable nuclear matrix elements by including 
the above refinements in the nucleon current in conjunction with the recent improvements 
of QRPA (renormalization effects due to Pauli principle corrections [^,^). In particular 
to see what effects, if any, the weak-magnetism and pseudoscalar coupling terms will have 
on the neutrino mass mechanism as well as on the Majoron emission mechanism of Oi>f3f3- 
decay. To this end we have performed calculations, which cover most of the nuclear targets 
of experimental interest C^Ge, ^^Se, ^^Zr, ^^^Mo, ^^^Cd, ^^^Te, ^^^Te, ^^^Xe, ^^^Nd). 

The paper is organized as follows: In Sect. II, the basic elements of the theory of the 
Oi//3/3-decay relevant to this work are presented. In the Sect. Ill the contributions coming 
from these induced currents to the Oz//?/5-decay amplitude are analyzed. Section IV summa- 
rizes the basic ingredients of the proton-neutron RQRPA (pn-RQRPA) method, which will 
be used for nuclear structure studies of the Oi//3/3-decay transitions. In Sect. V we discuss 
the calculation of nuclear matrix elements and deduce limits on lepton-number violating 
parameters. Finally, in Sect. VI we summarize the results and draw some conclusions. 



II. THEORY 
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A. Majorana neutrino mass mechanism 



We shall consider the Oi//3/3-decay process assuming that the effective beta decay Hamil- 
tonian acquires the form: 

= ^[ei,{l-i,)v,L]Ji^ + h.c., (1) 

where e and VeL are field operators representing electron and left handed electron neutrino, 
respectively. We suppose that neutrino mixing does take place according to 

^eL= E U^kXkL+ E U^kNkL, (2) 
k=light k=heavy 

where, Xfe (^fc) are fields of light (heavy) Majorana neutrinos with masses irtk (m^ << 1 

MeV) and (M^. >> 1 GeV), respectively, and f/^ is a unitary mixing matrix. In the 

first and second terms on the r.h.s. of Eq. (|]) the summation is only over light and heavy 

neutrinos, respectively. The fields Xk and Nk satisfy the Majorana condition: Xkik = C x!, 
— J' 

Nk^k = C Nf^ , where C denotes the charge conjugation and ^, ^ are phase factors. 

We assume both outgoing electrons to be in the Si/2 state and consider only Of Oj 
transitions are allowed. For the ground state transition restricting ourselves to the mass 
mechanism one obtains for the Oi//3/3-decay inverse half-life p|-|^,|6| JT3|JT^ , 



[r,%]-^ = Goil^^Ml^tV + Vr.Mi-y\'. (3) 

The lepton-number non-conserving parameters, i.e. the effective neutrino mass < > and 
rjj^ in Eq. (H) are given as follows: 

light heavy 

< > = E (U^'kf rrik, V. = {U^k? L ^, (4) 

k k ^ 

with rrip {rrie) being the proton (electron) mass. Gqi is the integrated kinematical factor 
PJ^. The derivation of the nuclear matrix elements associated with the exchange of light 
{M^<tl>) and heavy {M^^^^^y) Maj orana neutrinos is outlined in the next section. However, 
Eq. (P) apphes to any intermediate particle. 
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B. Majoron mechanism 



If the global symmetry associated with lepton number conservation is broken sponta- 
neously, the models imply the existence of a physical Nambu-Goldstone boson, called Ma- 
joron which couples to neutrinos: 

C<f,uv = J2 '^il5J^j (f)) Pij, Pij = J2 U,^ Ufp gap, (5) 

i<j a,f3=e,iJ,,T 

Here, Vi denotes both light Xi heavy Ni Majorana neutrinos. We remind that in analogy 
with Eq. (^) there is a unitary transformation for the right handed electron neutrino to the 
mass eigenstates Xk and Nk- 

yeR= E ^IXkR^ E Vek^^R. (6) 
k=light k=heavy 

where XkR = PRXk and NkR = PffiVfe (Pr,l = 1/2(1 ± 75)). 

Eq. d^) leads to Majoron production in the Oi//3/3-decay (Oi//5/30-decay) [Pt^|l3ip3| . We 
are interested in the light neutrino coupling and notice that the couplings are small in 
GUT models where the singlet neutrino is superheavy. We restrict our consideration of the 
Oi//3/50-decay only to light neutrinos ( rriij <^ q ^ pp ^ O{100MeV) ). Then the inverse 
half-life of the Oi//3/30-decay can be written 

[Tt;'n-' = \<9>\'m'j:!>\'GB. (?) 

Here < g > is the effective Majoron coupling constant 

light 

<9>=i: U^iU^^P^r (8) 
The explicit form of the kinematical factor Gb can be found in Ref. 0. 

III. THE EFFECTIVE TRANSITION OPERATOR 

Within the impulse approximation the nuclear current in Eq. (|T]) expressed with 
nucleon fields ^ takes the form 

6 



2m, 



-lu - 9A{q^)l^lb + 9p 



(9) 



where = {p — p')^ is the momentum transferred from hadrons to leptons {p and p' are 
four momenta of neutron and proton, respectively) and cr^" = (i/2)[7'^, 7*^]. gv{q^)i QMio^), 
gA{(f) and gp{q^) are real functions of a Lorenz scalar q^. The values of these form factors 
in the zero-momentum transfer limit are known as the vector, weak-magnetism, axial- vector 
and induced pseudoscalar coupling constants, respectively. 



A. The effective transition operator in momentum space. 

For nuclear structure calculations it is necessary to reduce the nucleon current to the 
non-relativistic form. We shall neglect small energy transfers between nucleons in the non- 
relativistic expansion. Then the form of the nucleon current coincides with those in the 



Breit frame and we arrive at 124 



J'{x) = j2r+[g^'j'{q^)+g^'Jt{q^mS-r^), k = 1,2,3, (10) 

n=l 



with 



J\q')=gv{q'), Jn{q')=gM{qY-^+9A{q')a-gp{q')^-^. (11) 

r„ is the coordinate of the rath nucleon. 

For the form factors we shall use the following parameterization: gviq"^) = 
gv/i^ + q^/K?, gAiiq^) = if^p- f^n)gviq^), gA{q^) = gA/{^ + q"^ and the induced 
pseudoscalar coupling is given by the partially conserved axial-vector current hypothesis 
(PCAC) [H 

2 

777 

gp{q') = 2m,gA{q')/{q' + m^)(l - ), (12) 

^^A 



where gv = I, gA = 1-254, (/ip-/i„) = 3.70, = 0.71 {GeVY [gUf and A^ = 1.09 GeV ^ 
In previous calculations only one general cut-off Ay = A^ ^ 0.85 GeV was used. In this 
work we take the empirical value of A^ deduced from the antineutrino quasielastic reaction 



Vfj_p — > ji'^n. A larger value of the cut-off is expected to increase slightly the values of 
the corresponding nuclear matrix elements. It is worth noting that with these modifications 
of the nuclear current one gets a new contribution in the neutrino mass mechanism, namely 
the tensor contribution. 

As we have already mentioned in the introduction, momentum dependent terms, in 
particular the weak-magnetism term, have been considered previously in the (3(3 decay by 
Tomoda et al. and Pantis et al [§, but in connection with the 77-term. This term is 



proportional to the mixing between the vector bosons Wl and Wr-, which mediate the left 
and right handed weak interaction, respectively. They are dominant since, due to their 
momentum structure, they can proceed via the s-wave electron wave function, while the 
standard terms in this case require p-wave electron wave functions. The pseudoscalar term 
is not accompanied by parity change and thus it is not important in the extraction of 77. 
To our knowledge this term has not been considered in connection with the usual light 
Majorana neutrino mass term of the Oz//3/3-decay. 

Under the PCAC hypothesis [see Eq. (|T2|)] the two body effective transition operator 
takes in momentum space the form 

^ = r+T+{-hF + hcT (7i2 - hTSi2) (13) 

where the three terms correspond to Fermi (F) , Gamow- Teller (GT) and Tensor (T). One 
finds that 

5'i2 = 3(cti • qa2 ■ q) - (T12, (T12 = ^1 ■ ^2- (14) 
hp = 9v{q^) 

.2^ 2,^2^ n 2 q' , 1 , q^ 1 , 2 gUq')q' 



hcAV) = 9iir) [1-77 + o i^T^r ] + 



3 g ^ + ml 3 q'^ + ml 3 Am^ 

The exact results will depend on the details of the nuclear model, since the new operators 
have different momentum (radial) dependence than the traditional ones and the tensor 
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component is entirely new. We can get a crude idea of what is happening by taking the 
above average momentum (g)=100 MeV. Then we find that the GT matrix element is 
reduced by 22%. Then assuming that T matrix element is about half the GT one, we find 
that the total reduction is 28%. This will be compared below with the results of our detailed 
calculations. 



B. The effective transition operator in coordinate space. 

The nuclear matrix elements entering the half-hfe formula of 0i//9/9-decay process take 
now the form: 

Km.>,v^ = M(rv + M'mm + Mij, + Mip + M'pp (16) 

with / = light, heavy. The partial nuclear matrix elements M^y, Mljj^, M^^, Mpp and 
Mj^p have their origin from the vector, the weak-magnetism, the axial, the pseudoscalar 
coupling and the interference of the axial-vector and pseudoscalar coupling, respectively. 
They can be expressed in relative coordinates by using second quantization. We end up 
with formula 

Kpe =< Kpe-Firu) + Hlpe-GTiri2)ai2 + i^4pe-r(n2)^12) > (17) 

with type = VV, MM, AA, PP, AP and 

_ _ I I - _ "^12 

ri2 — Ti — r2, ri2 — i'i2 , t^i2 — — , 

5'i2 = 3(cti • fi2)(a2 • fi2) - (712, (712 = CTi • (72- (18) 

ri and r2 are coordinates of the beta decaying nucleons. The form of the matrix element 
< 0{1, 2) > within the pn-RQRPA will be presented in the next section. 

The light and heavy neutrino-exchange potentials H^typ^-K^^^ {^'12) {K = F,GT,T) have 
the following forms, 

1 2 R 



rrUght ( 2 R smjqru) 

type-KK 12) ^^2 ]^ q + Ef - (El,, + EU)/2 

1 2 R f°° 

Htype-Kiri2) = o — / sin(gri2)/it,pe-i^ (g') Q dq (20) 

m„mp, 77 g A rio jo 



Here, Eg ^ , ^ and EJ' are respectively the energies of the initial, final and intermediate 
nuclear states. R = tqA^^^ is the mean nuclear radius, with tq = 1.1 fm. The relevant 
couplings are: 

hvv{q') = -9UQ') 

) - 3 Ami ' ^ - 3 4m2 ' 

) = 3 ' ^^^-^^^ ^ = "3 4m^ ' 

(-^2^ '2 gA{q^)gp{q^)q^ . 2^ QAjq ^)gp{q ^)q ^ 

) = "3 2^, ' ^^^-^^^ ^ = 3 2^;^; ■ ^^^^ 

The tensor form factor includes a sign change going from momentum to coordinate space. 

The full matrix element is of the form: 

T\Aight 

M'^tl> = + Mj^f + M^''''- (22) 

9a 

We see that the Fermi component is unchanged, the Gamow- Teller is modified and the tensor 
component appeared due to the new terms. 





-GT 


hAA- 


-GT 


hpp- 


-GT 


Hap- 


-GT 



IV. NUCLEAR STRUCTURE INGREDIENTS 

As we have mentioned above, we would like to evaluate the changes in the Oz//?/3-decay 
nuclear matrix elements due to the modifications of the nuclear current introduced above, 
relevant for the neutrino mass mechanism. It is clear that the Oz//3/3-decay is a second 
order process in the weak interaction and, thus, the corresponding nuclear matrix elements 
require the summation over a complete set of intermediate nuclear states. Even though 
the construction of these states is not needed, a closure approximation with a reasonable 
average energy denominator is very accurate [p^-|28[], the initial and final states of the nuclear 
systems, which can undergo double beta decay, are not easy to construct, since these nuclei 
are far removed from closed shells. Thus the introduction of additional approximations is 
necessary. 

10 



Thus at this point we will reduce the computational difficulty in evaluating the effects 
of the above mentioned modifications, by using the proton-neutron Quasiparticle Random 
Phase Approximation (pn-QRPA) pO|-p^], which is an approximation to solve the nuclear 
many-body problem. Admittedly the intermediate states must be explicitly constructed in 
this case, but in this scheme it is a a simple matter to include even a large number of such 
states, if necessary. 

The crucial simplifying point of the QRPA is the quasiboson approximation (QBA) 
assuming nuclear excitations to be harmonic. It leads to a violation of the Pauli exclu- 
sion principle. The drawback of this approximation is that with increasing strength of the 
nucleon-nucleon force in the particle-particle channel the QRPA overestimates the ground 
state correlations and the QRPA solution collapses pO| , pT| , p4| , |35| . 



The renormalized QRPA (RQRPA) p0|j21[| overcomes this difficulty by taking into ac- 
count the Pauli exclusion principle in a more proper way by using the renormalized QBA. 
In this paper we calculate the nuclear matrix elements of the Oz//5/5-decay within the proton- 
neutron RQRPA (pn- RQRPA) [^,^, which is an extension of the pn-QRPA by incorpo- 
rating the Pauli exclusion principle for fermion pairs. For studying the relative importance 
of the new induced terms, which is the main thrust of our paper, the inclusion of other 
refinements like p-n pairing, which is computationally very involved, is not essential. 

Furthermore p-n pairing is rigorously incorporated in the BCS ansatz only for the T=l 
states, while T=0 p-n pairing effects are implicitly taken into account. This procedure seems 
to avoid the collapse of the QRPA within the physical region of the Hamiltonian but we are 
not sure whether in some cases it does not produce more ground state correlations which can 
lead to strong cancellations in the matrix element. This might have been the case of ^°°Mo 
in our earlier work H] in which one should have expected a minor influence of p-n pairing. In 



fact if the collapse of the QRPA reflects a nearby phase transition ||3^, i.e., a change of the 
ground state from being dominated by T=l pair correlations to being dominated by T=0 
pair correlations, further work needs to been done to be sure about the p-n pairing. On the 
other hand the renormalized solution does not lead to a collapse of the QRPA for physical 
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values of the proton-neutron interaction strength and is tractable from the computational 
point of view. 

The pn-RQRPA excited states |m, JM > are of the form [pO| , pl| 



|J-Mm) = Q7,V|0W 

= E[X(;n,j^)A\pn, JM) + F(^„,j.)i(pn, JM)]|0+p^), (23) 

pn 

where XJ^^ and j^-j are free variational amplitudes, respectively, and 

A\pn, JM) = Y: C^i^,,„rn„4-,«L„, Mpn, JM) = (-l)^-^A(pn, J - M). (24) 

mp,m„ 

Here, a+^^ (armr , t = P,n) is the quasiparticle creation (annihilation) operator for spherical 
shell model states, which is related to the particle creation and annihilation (c:^^^ and c^^^ , 
T = p,n) operators by the Bogoliubov-Valatin transformation: 



(25) 



where the tilde ~ indicates time reversal {drm-r = (~l)"'^~™'^'^r-mr)- '^^^ label r designates 
quantum numbers n^, It-^Jt- The occupation amplitudes u and v and the single quasiparticle 
energies are obtained by solving the BCS equation. 

In the pn-RQRPA the commutator of two-bifermion operators fulfill the following relation 
(renormalized QBA) 

{Q-tiPA\[Mpn.JM),A^{p'n\JM)]\Qlp^)=5,p,p,5nn' x 

(l - — < Ojp^|[a+ap]oo|Ojp^ > - — < 0^p^| [a+a„]oo|0^p^ > |, 

^ Jl Jk ^ 

V ' 

T-^prijJ^ 

(26) 



with jp = y2jp + 1. If we replace |0^p^ > in Eq. ( p6D with the uncorrelated BCS ground 
state, we obtain the quasiboson approximation (i.e. Vpnj^ = 1), which assumes that pairs 
of quasiparticles obey the commutation relations of bosons. We note that it is convenient 
to introduce amplitudes 
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Y — T) ' Y"^ V — T) ' V"^ (07\ 

which are orthonormahzed in the usual way 

Smm' = X! i-^ {pn,J^)^ (pn,J^) ~ ^ (pn,J^)^ (pn,J^)) ■ (28) 



pn 

One can show that 

< J'^rrii II [c+Cn]j II 0+ > 
V2J + 1 



< 0/ II [c+5n]j II J^'ruf > 



[Up ^{pn,J^) ^ "p "n {pn,J^))^ ^pn,J-^i ) 



[Vp ^ {pn,J-^) + ^p ^ {pn,J^))\l ^pn,J^- 



V27TT 

The index i (f) indicates that the quasiparticles and the excited states of the nucleus are 
defined with respect to the initial (final) nuclear ground state |0j^ > (|0j >). The forward- 
and backward- amplitudes ^^^^j^) and jtt) and the energies of the excited states fi^'? = 
E^T^ — Eg s are obtained by solving the non-linear set of RQRPA equations for the initial 
nucleus (A,Z) [^,^. By performing the RQRPA diagonalization for the final nucleus 
(A,Z=2) we obtain the amplitudes X^^^^j-^-j and j^r) and the eigenenergies ^I'JI = EJ'J — 
El^ of the RQRPA state \J^mf). 

Within the pn-RQRPA the Oz//?/?-decay matrix elements given in Eqs. (|T^ take the 
following form: 

jp in J 
in' jp' J 

<p{l),p\2)-J\f{n2)rtT^Ol^,{l2)f{n2)\n{l),n\2)-J>^ 
< 0/ II [cpn']j II rnif >< rmfirnii >< rnii II [c+c„]j II 0+ > . (31) 

Here, Olyp^{12) represents the coordinate and spin dependent part of the two body transition 
operators of the Oz//5/5-decay nuclear matrix elements in Eq. ( ^I]) 



MLe= E (-)^"+V+^+^(2J+l) 



type 

pnp' 



X 



0\i\2) = H'in,) + H'^{n,)a,2 + H'{ru)S,2. (32) 
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The short-range correlations between the two interacting protons and p'(2)) and neu- 

trons {n{l) and n'{2)) are taken into account by the correlation function f{ri2) in the 
non-antisymmetrized two-body matrix element in Eq. (pi]), firu) is given as follows: 

/(na) = 1 - e-"'i2(l - brl^) with a = 1.1 fm^ and b = 0.68 fm^ (33) 

For the overlap matrix of intermediate nuclear states generated from the initial and final 
ground states we write ||35|| : 



pn 



(««4^) + vl:^vi/'>){n^n^^ + t;«4^)). (34) 



V. CALCULATION, DISCUSSION AND OUTLOOK 



In order to test he importance of the new momentum dependent terms in the nucleon 
current, we applied the pn-RQRPA to calculate the Oi//?/3-decay of the A = 76, 82, 96, 100, 
116, 128, 130, 136 and 150 systems. To this end the considered single-particle model spaces 
both for protons and neutrons have been as follows: i) For A=76, 82 the model space consists 
of the full 2 — Ahuj major oscillator shells, ii) For A=96, 100, 116 we added to the previous 
model space I/5/2, I/7/2, O/19/2 and 0/iii/2 levels, iii) For A=128, 130, 136 the model space 
comprises the full 2 — 5huj major shells, iv) For A=150 the model space extends over the 
full 2 — 5hu shells plus the Oiii/2 and 0213/2 levels. 

The single particle energies were obtained by using a Coulomb corrected Woods Saxon 
potential. The interaction employed was the Brueckner G-matrix which is a solution of the 
Bethe-Goldstone equation with the Bonn one-boson exchange potential. Since the model 
space considered is finite, the pairing interactions have been adjusted to fit the empirical 



pairing gaps according to ^Tj . In addition, we renormalize the particle-particle and particle- 
hole channels of the G-matrix interaction of the nuclear Hamiltonian H by introducing the 
parameters Qpp and Qph, respectively. The nuclear matrix elements listed in Tables | and H 
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have been obtained for gph = 0.8 and gpp = 1.0. With respect to the gpp we wish to make 
the following statement: Our numerical results do not show significant variations (do not 
exceed 20 %) in the physical region of gpp (0.8 < gpp < 1.2). 

A detailed study of Fermi, Gamow- Teller and Tensor contribution to the full nuclear 
matrix element M<^^„*> in Eq. (E^f) for the two representative Oi//5/?-decay nuclei ^^Ge and 



^^^Te is presented in Table |. One notices significant additional contributions to GT (AA 
and PP) and tensor (AA and PP) nuclear matrix elements coming from the induced current 
terms. By glancing at the Table | we also see that Fermi and GT matrix elements are 
strongly suppressed due to the nucleon short-range correlations. 

The relative importance of the different contributions to the nuclear matrix elements 
M<^''^*> and M^^^'^y [see Eq. (||)] is displayed in Fig. |l] for the A=76 and 130 nuclear 
systems. As expected from our discussion in Section III we see that for light neutrinos the 
pseudoscalar term, in particular the AP contribution, is important. It is in fact as important 
as the usual vector contribution but in the opposite direction (see Table fT|) . Our calculations 
verify our above estimate, i.e. the new terms in the hadronic current, and in particular the 
induced pseudoscalar term, tend to increase the average neutrino mass | < > \ and 
average Majoron coupling < g > from experiment by about 30 %. They are much more 
important in the exchange of heavy neutrinos leading the suppression of M^*^"^^ by about 
factor of 3-6 (see Table H). The contributions from previously neglected M^^^ and M^p^'^ 
to M^'^"''"y are much more important as that from My^"^. A large value of M^^^, which 
has its origin in weak-magnetism, indicates that the average neutrino momentum is large, 
i.e. about the order of magnitude of the nucleon mass. 

We present in Fig. ^ the nuclear matrix elements M<mJ> and M^*^"^^ calculated within 
pn-RQRPA for the A = 76, 82, 96, 100, 116, 128, 130, 136 and 150 nuclear systems. We 
see that the inclusion of the induced pseudoscalar interaction and of weak-magnetism in the 
calculation results in considerably smaller nuclear matrix elements for all nuclear systems. 
The numerical values of M<m^*> and M^'^"-'"^ can be found in Table ||. The largest matrix 
elements of M<m^*> for A=150, 100 and 76 are 3.33, 3.21 and 2.80, respectively. For A=150, 
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76 and 82 the largest values of M^^'^^?' are 35.6, 32.6 and 30.0, respectively. We notice that the 
A=136 system has the smallest nuclear matrix elements: M<^^*> = 0.66 and M^*^™^ = 14.1. 
We suppose that it is connected with fact that ^^^Xe is a closed shell nucleus for neutrons 
(N=82). The sharp Fermi level for neutrons yield smaller Oz//3/3-decay matrix elements. We 
note that in calculating the matrix elements involving the exchange of heavy neutrinos, the 
treatment of the short-range repulsion and nucleon finite size is crucial. We have found 
that the consideration of short range correlation effects reduces the values of M^*^™^ by 
about factor of 20-30. As we mentioned already the nucleon finite size has been taken into 
account through the phenomenological formfactors and the PCAC hypothesis. However, the 
choice of the formfactor can influence the results significantly as it was manifested in Ref. 
T8| performing the calculations with both phenomenological and quark confinement model 



formfactors. 

The limits deduced for lepton-number violating parameters depend on the values of 
nuclear matrix element, of the kinematical factor and of the current experimental limit for 
a given isotope [see Eqs. (^) and (0)]. It is expected that the experimental constraints on 
the half-life of the 0z//3/3-decay are expected to be more stringent in future. Thus there is 
useful to introduce sensitivity parameters for a given isotope to the effective light and heavy 
Majorana neutrino mass and Majoron signals, which depend only on the characteristics of 
a given nuclear system. There are as follows: 



< m,, > , 10" 



C<m.>iY) = 10^ |M^t*J V^Goi year, — < , 



10^4 years 



C (Y) = 10« \M^Zl I V^oi year, r^, < 



10 



-6 



10^4 years 



C<,>(r) = 10« iM^f^Vl ^Gb year, < g >< 



10-^ 



rpOu—exp ' 
^1/2 



10^4 years 

rpOu—exp 
-^1/2 



(35) 



C<.>\ 

The normalization of 10^^ years was chosen so that the C's are of order unity. The numerical 
values of these parameters for the nuclear systems considered in this work are listed in Table 
|T| and can be used in predicting desired limits on lepton number violation with changing 
experimental data (see the above expressions). These characteristics estimate also prospects 
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for searches of hght and heavy Majorana neutrinos and of the Majoron. The larger values 
of these parameters determine those Oi//5/5-decay isotopes, which are the most promising 
candidates for searching for the corresponding lepton number violating signal. They give 
information on the requirements on a Oi//5/5-detector. 

Fig.0 shows that the most sensitive isotope to all three studied lepton number violating 
parameters is ^^^Nd. It is mostly due to the large phase space integral and to some extent 
also to the larger nuclear matrix elements. We caution the reader here that our matrix 
elements are not model independent in the sense that all the nuclei considered in this work 
are treated as spherical . The nucleus ^^^Nd, is deformed and our results may not be the 
same quantitatively, were we to perform calculations taking into account effects of nuclear 
deformation. 

The purpose of this work, however, is to study the effects of the induced currents. It is 
thus reasonable to do this by comparing calculations within the same model. This means 
that the model itself will play a minor role, if any, in the investigation of such effects. We are 
thus satisfied that for light neutrinos the effect is almost the same throughout the periodic 
table. In addition, as we see from Fig. the inclusion of weak-magnetism and induced 
pseudoscalar coupling in our calculations leads to a significant reduction of the parameters 
of C<™.>(n,C,,(n and (<,>{¥). 

The present experimental situation in terms of the accessible half-life and the correspond- 



ing upper limit on < mj, >, i]^ and < > is given in Table |T|. Thus, the most restrictive 
limits follows: 



< >'^'' < 0.62 eV, pCe, Ref. 38] 

< g < 6.9 X lO-l [^2^Te, Ref. 43]. (36) 

The sensitivity of different experiments to < m^, >, < rj^ > and < g > is drawn in 



Fig. ^. Currently, the Heidelberg-Moscow experiment offers the most stringent limit 
for effective light and heavy Majorana neutrino mass and the ^^^Te experiment for 
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effective Majoron coupling constant. By assuming < m,y >=< ^best^ _ ^hest ^.^^ 



N 



< g >=< g >best gqs. (D and (g) we calculated half-lifes of the Oz//3/5-decay T^J^'^^iK 
nil, >^^^^)^ Ty2~'^'^ iv'^r^'^^) Ty2~'^'^'^{< 9 >^'^^*) for nuclear systems of interest using specific 
mechanisms with the " best " parameters. These corresponding numerical values are listed 
in Table ^ and shown by open bars in Fig. ^. Since the quantities < rrii, >, 1]^^, < g > 
depend only on particle theory parameters these quantities indicate the experimental half- 
life limit for a given isotope, which the relevant experiments should reach in order to extract 
the best present bound on the corresponding lepton number violating parameter from their 
data. Some of them have a long way to go to reach the ^^Ge target limit. 

At present most attention is paid to the light Majorana neutrino mass because of the 
experimental indications for oscillations of solar (Homestake Kamiokande [Q, Gallex 
)1[ and SAGE [0), atmospheric (Kamiokande |5^, 1MB and Soudan [0], Super- 



Kamiokande experiments [Q) and terrestrial neutrinos (LSND experiment [0). One can 
use the constraints imposed by the results of neutrino oscillation experiments on < rriu >. 
The predictions differ from each other due the different input and structure of the neutrino 
mixing matrix and in particular the assumed Majorana condition phases. Bilenky et al 
have shown that in a general scheme with three light Majorana neutrinos and mass hierarchy 
I < rriu > I is smaller than 10~^ eV. In another study outlined in Ref. |^9| the authors end 



up with I < rrii, > \ ^0.14 eV. Bednyakov, Faessler and Kovalenko considered neutrino 
oscillations within the minimal supersymmetric standard model with R-parity breaking. 
They showed that Super-Kamiokande atmospheric data are compatible with | < mj, > | < 
0.8 X 10~^ eV [|^. One sees that, the current limit on < > in (p6D is quite a bit higher 



than the neutrino oscillation data. 

There is a new experimental proposal for measurement of the Oz//5/5-decay of ^^Ge, which 
intents to use 1 ton (in an extended version 10 tons) of enriched '^^Ge and to reach the 
half-life limit Ty^^""" > 5.8 x 10^^ and T°/7'^'^ > 6.4 x lO^^ after one ( and 10 years) of 
measurements, respectively. From these half-life values one can deduce [see Eq. and 
Table |l| the possible future limits on the effective light neutrino mass 2.7 x 10~^ eV and 



8.1 X 10~^ eV, respectively. 

By comparing the above limits with those advocated by the neutrino oscillation phe- 
nomenology we conclude that the GENIUS experiment will be able to reach similar limits, 
provided, of course, that the neutrinos are Majorana particles. We remind that there is also 
a plethora of other Oi//?/?-decay mechanisms predicted by GUT's and SUSY. One can show, 
however, that their presence implies that the neutrinos are massive Majorana particles even 



if the mass mechanism is not dominant l6^j6^ . Certainly, the experimental detection of the 
Oi//3/?-decay process would be a major achievement with important implications on the field 
of particle and nuclear physics as well as on cosmology. 

VI. CONCLUSIONS 

The contributions coming from the induced currents at the nucleon level, such as the 
weak- magnetism and induced pseudoscalar coupling on the mass mechanism for the Qv(3(3- 
decay transitions has been studied. The needed nuclear matrix elements, associated with 
the light and heavy Majorana neutrinos as well as the Majoron emission mechanisms, have 
been obtained in the context of pn-RQRPA, which is known to produce results more reliable 
than the standard QRPA. Our results are shown in Figs. and listed in Tables |l-|l|. One 
can see that the modification of the nuclear current due to the weak-magnetism and induced 
pseudoscalar coupling is important and results in considerable reductions of the Oz//5/5-decay 
matrix elements. For the light neutrino exchange this reduction amounts to about 20 — 30% 
for all nuclei considered. The reductions for the heavy neutrino exchange are even more 
significant with factors ranging from 4 to 6. 

The derived upper limits on < niy >, rjj^ and < g > from the current experimental limits 
of the Oi//3/3-decay lifetime for A = 76, 82, 96, 100, 116, 128, 130, 136 and 150 are listed in 
Table |T1|. This makes the extracted limits of the lepton number violating parameters less 
stringent yielding < ^best^ q q2 ^ ^best < ^ q ^ < g ^best < g g ^ iq-s 



deduced from the ^^Ge PB[ and ^^^Te P3] data. Further, we introduced and evaluated useful 
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sensitivity parameters for various lepton number violating signals for some nuclei of interest, 
which might be helpful in planning future Oz//3/5-decay experiments. 

The value of 1]^ extracted is, of course, associated with heavy Majorana neutrino. It 
can, however, be applicable in other processes involving the exchange of heavy particles, 
provided that the momentum structure of the relevant operators is not very different from 
that in Eqs. (|19D-(P). 



Admittedly there is a rather large spread between the calculated values of nuclear matrix 



elements within different nuclear theories [see e.g. recent review articles [|T1,Q], which could 
be considered as a measure of the theoretical uncertainty. Between some of them there is no 
objective way to judge which calculation is correct. However, one can argue that the RQRPA 
method offers more reliable results than the QRPA primarily because of the collapse of the 
QRPA solution and the strong sensitivity of the QRPA results to the strength of particle- 
particle force. The only advantage of the QRPA over RQRPA is that it fulfills the Ikeda 
sum rule. However, the meaning of this fact is questionable because it is so close to the 
collapse of the QRPA, where the obtained solution is far from being realistic. In the present 
calculations we are using the pn-RQRPA and we take into account also additional nucleon 
currents effects. Thus we consider the results of this paper more reliable in respect to pn- 
QRPA results of our and other groups. In this work we did not deal with the problem of 
the proton-neutron pairing. The effects of proton-neutron pairing within the renormalized 



QRPA have been discussed for some nuclei of interest in Refs. ^4 . 

Be that as it may, we find that in the case of light neutrino the momentum dependent 
terms in the nucleon current cause a more or less uniform reduction of the nuclear matrix 
elements by approximately 30% throughout the periodic table. We expect a similar reduction 
in almost any nuclear model. This will cause a corresponding increase of the extracted values 
for the neutrino mass. 

We thus conclude that, with the best nuclear physics input, the extracted average neu- 
trino mass is low, but quite a bit higher than that deduced from the present neutrino 
oscillation experiments. It will reach, however, similar levels, if the planned experiments 
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reach the sensitivity aimed at by the GENIUS experiment [^. In any case the neutrino os- 
cillation data can neither set the absolute scale of the mass nor decide whether the neutrino 
is a Majorana particle. The latter issue can be decided only by the Oi//3/5-decay. 
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TABLES 

TABLE I. The Fermi, Gamow-Teller and Tensor nuclear matrix elements for the light Majo- 
rana neutrino exchange of the Oi//3/3-decay of ''^Ge and ^^^Te with and without consideration of 
short-range correlations (s.r.c). 

transition s.r.c. Gamow-Teller Tensor Mp^* M^f^* M^'^''* 





AA 


AP 


PP 


AP 


PP 








76(5g ^76 without 


5.132 


-1.392 


0.302 


-0.243 


0.054 


-2.059 


4.042 


-0.188 


with 


2.797 


-0.790 


0.176 


-0.246 


0.055 


-1.261 


2.183 


-0.190 


isorp^ ^130 xe without 


4.158 


-1.173 


0.258 


-0.329 


0.074 


-1.837 


3.243 


-0.255 


with 


1.841 


-0.578 


0.131 


-0.333 


0.075 


-1.033 


1.397 


-0.258 
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TABLE II. Nuclear matrix elements for the hght and heavy Majorana neutrino exchange modes 
of the Oi^/3/3-decay for the nuclei studied in this work calculated within the renormalized pn-QRPA. 

Goi and Gb are the integrated kinematical factors for the 0+ 0+ transition. C<m^>(^)) CtinO^) 
and C<g> (X) denote according to Eq. (34) the sensitivity of a given nucleus Y to the light neutrino 

mass, hcavj- uculriuo mass and Majorou signals, rosi")eclively. 







(/3/3)oi/ - decay : 0+ - 


0"*" transition 








M P 

IVl. JZj. 














J. c- 


136 Xp 


150 AT J 
IM U 








light Maj 


orana neutrino (1= 


dight) 










0.80 


0.74 


0.45 


0.82 


0.50 


0.75 


0.66 


0.32 


1.14 




z.ou 


2.66 


1.54 


3.30 


2.08 


2.21 


J. .o^ 


U. I u 


o.o / 


IVlpp 


U./o 


0.22 


0.15 


0.26 


0.15 


0.24 




nil 

U.ii 


n Qc; 






-0.98 


-0.65 


-1.17 


-0.69 


-1.04 


n Qi 


n /IS 


-i.OO 


Myy + Mj^j^ 


Q fin 


3.40 


1.99 


4.12 


2.58 


2.96 


z.ou 


1 no 


A t^l 
4.01 




2.80 


2.64 


1.49 


3.21 


2.05 


2.17 


1.80 


0.66 


3.33 






heavy Majorana neutrino (1= 


heavy) 








]\ /fl 


23. y 


22.0 


16.1 


28.3 


17.2 


25.8 




ion 

13.9 


39.4 




-55.4 


-51.6 


-38.1 


-67.3 


-39.8 


-60.4 


-54.5 


-31.3 


-92.0 


Ka 


106. 


98.3 


68.4 


123. 


74.0 


111. 


100. 


58.3 


167. 


Ml,p 


13.0 


12.0 


9.3 


16.1 


9.1 


14.9 


13.6 


7.9 


23.0 


Kp 


-55.1 


-50.7 


-41.1 


-70.1 


-39.0 


-64.9 


-59.4 


-34.8 


-101. 


Kv + Ka 


130. 


120. 


84.5 


151. 


91.1 


137. 


123. 


72.3 


206. 


In 


32.6 


30.0 


14.7 


29.7 


21.5 


26.6 


23.1 


14.1 


35.6 



sensitivity to neutrino mass signal 
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Goi X lO^^y 



7.93 


35.2 


73.6 


57.3 


62.3 


2.21 


55.4 


59.1 


269. 


2.49 


4.95 


4.04 


7.69 


5.11 


1.02 


4.24 


1.60 


17.3 


2.90 


5.64 


3.98 


7.10 


5.36 


1.25 


5.45 


3.43 


18.5 



sensitivity to Majoron signal 
Gb X 7.40 62.3 159. 106. 104. 0.59 79.6 82.8 640. 

C<g>(y) 2.41 6.59 5.93 10.5 6.60 0.53 5.08 1.90 26.7 
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TABLE III. The present state of the Majorana neutrino mass and Majoron searches in /3/3-decay 
experiments. T^J|'~'''^ (present) and T-[^J|'~^'^''^ (present) are the best presently available lower limit 
on the half-life of the Oi^/?/3-decay and Oz^/?/3i;^-decay for a given isotope, respectively. The corre- 
sponding upper limits on lepton number non-conserving parameters < rrii, >, < g > and rj^ are 
presented. T^J^-^^iK >^^^*), rf/r°''(^^''*) and Tl'J^-'^'"^ {< g >^^^*) are calculated half-lifes of 
Oz^/3/3-decay assuming < rrii, >=< rriy >^^'^*, = r/^^''* and < g >=< g >''^''*, respectively. Here, 



< my >^est^ Q g2 eV 7?^^''* = 1.0 x 10"^ and < g >''e^«= 6.9 x 10"'' are the best limits deduced 



J)est 



best_ 



from the ™Ge [38] and ^^^Te [43] experiments. 



128r 



Nucleus 



76 



Ge 



82 



Se 



96 



Zr 



100 



Mo 



116 



C7d 



''(present) [y] 

Ref. 

<mu> [eV] 



Tl^r"'{< my >^-*) [y] 



1.1 X 10 



25 



0.62 

1.1 X 10^^ 

1.0 X 10"^ 

1.1 X 10^5 



2.7 X 10^2 3.9 X 10^^ 5.2 x 10^^ 2.9 x 10^2 



6.3 

2.8 X 10^4 
1.1 X 10-6 

2.9 X lO^'' 



203. 

4.2 X 10^4 
4.0 X 10-5 
5.8 X lO^'' 



ill 
2.9 

1.2 X 10^4 
6.2 X lO-'^ 



1.8 X 10 
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5.9 

2.6 X 10^4 
1.1 X 10-6 



3.2 X 10 



24 



j,exp 0'^'/' (^present) [y] 
Ref. 



<9> 

n^r'^^^rf^) [y] 



7.9 X 10^1 

m 

4.7 X lO"'^ 
3.7 X 10^3 



1.6 X 10^1 3.9 X 10^9 



@ 



3.8 X 10-4 

4.9 X 10^2 



2.7 X 10-3 
6.0 X 10^2 



5.4 X 1021 


1.3 X 10-4 

1.9 X 1022 



1.2 X 1021 

in 

4.4 X 10-4 
4.9 X 1022 



Nucleus 



128 



Te 



130 



Te 



136 



Xe 



150 



Nd 



"'^ (present) [y] 

Ref. 

<my> [eV] 
Tl^r^'iK my >^-*) [y] 



7.7 X 1024 8.2 X 102i 4.2 x 10 



23 



m 

1.8 

6.6 X 1025 
2.9 X lO-'^ 
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5.9 X 10^5 3.1 X W'^^ 7.9 x lO^'' 2.7 x 10^^ 



T^jr'-'Hwesent) [y] 
Ref. 

<g> 



7.7 X 2.7 X 10^1 1.4 X 10 
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3.8 X lO^'' 
8.2 X 10^2 



4.5 X lO"'' 
5.9 X 10^3 



ii 

2.2 X 10^-^ 
3.0 X 10^1 



30 



FIGURES 

FIG. 1. Calculated light and heavy neutrino exchange Oz//3/3-decay nuclear matrix elements 

for the A=76 and 128 systems. The partial matrix elements Myy-, Maa-. MmM: Mpp and Map 
originate from vector, axial-vector, weak magnetism, pseudoscalar coupling and the interference of 
the axial- vector and induced pseudoscalar coupling, respectively. M<^^> and M^^ are Oi//3/3-decay 
matrix elements associated with the < nii, > and r/j^ parameters, respectively. 

FIG. 2. Calculated nuclear matrix elements M<m^>, M^^^ [see Eqs. (3) and (13)-(18)], sen- 
sitivities C<mi,>5 Cvn C<ff> foi" the experimentally interesting A=76, 82, 96, 100, 116, 128, 
130, 136 and 150 nuclear systems. C<mi,>) Cv^ ^^'^ C<9> sensitivities to light neutrino mass, 
heavy neutrino mass and Majoron signal [see Eq. (34) ], respectively. The open and black bars 
correspond to results obtained without and with the inclusion of the pseudoscalar interaction and 
of weak-magnetism. 

FIG. 3. The sensitivity of different experiments to the lepton-number violating parameters 
< m,y >, r/^ and < g > are illustrated by histograms on the left side. The best presently available 
lower limits on the Oz//3/3-decay half-life ^'^ and Oz//3/30-decay T'^^^ ^'^'^ displayed by black 
bars in histograms on the right side. The open bars in these histograms indicate the half-life limits 

j^exp-ou^^ T;^/r°''«''*) and rf/f-°^'^(< g >^^^'*) to be reached by a given experiment 

to reach the presently best limit on < m,^ >, r]j^ and < g >, respectively. 
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